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Anyons exist as point like particles in two dimensions and carry braid statistics which enable 
interactions that are independent of the distance between the particles. Except for a relatively 
few number of models which are analytically tractable, much of the physics of anyons remain still 
unexplored. In this paper, we show how U(l)-symmetry can be combined with the previously 
proposed anyonic Matrix Product States to simulate ground states and dynamics of anyonic systems 
on a lattice at any rational particle number density. We provide proof of principle by studying 
itinerant anyons on a one dimensional chain where no natural notion of braiding arises and also on 
a two-leg ladder where the anyons hop between sites and possibly braid. We compare the result of 
the ground state energies of Fibonacci anyons against hardcore bosons and spinless fermions. In 
addition, we report the entanglement entropies of the ground states of interacting Fibonacci anyons 
on a fully filled two-leg ladder at different interaction strength, identifying gapped or gapless points 
in the parameter space. As an outlook, our approach can also prove useful in studying the time 
dynamics of a finite number of nonabelian anyons on a finite two-dimensional lattice. 


I. INTRODUCTION 

Anyons are point-like (quasi)particles which exist only 
in two-dimensional systems and have richer exchange 
statistics than bosons or fermions. One of the main in¬ 
terests in anyons is in their application to implement¬ 
ing fault-tolerant (topological) quantum computation. 
Anyons have also garnered a substantial theoretical inter¬ 
est since they are proposed to exist in systems as diverse 
as fractional quantum Hall systems and two-dimensional 
spin liquids,one dimensional nanowires,and 
ultra-cold atoms in optical lattices.Recent experiments 
showing evidence for Majorana edge modes (i.e. Ising 
anyons) in nanowires^° might bring us closer to work¬ 
ing with anyons in the lab, with far-reaching scientific 
and technological applications. 

One dimensional chains of static SU(2)fc anyons with a 
local antiferromagnetic Heisenberg-like interaction have 
been studied extensively since, for example, they are crit¬ 
ical and realize all minimal models of conformal field the¬ 
ories (CFTs).^^ It is also natural to ask whether interest¬ 
ing states and phases appear in anyon models where the 
anyons are allowed to hop on a lattice and braid around 
one another. Braiding pairs of anyons generally trans¬ 
forms the anyonic state in a non-trivial way, in contrast 
with bosons and fermions which merely pick up a fac¬ 
tor of ±1. For anyons, braiding is a topological interac¬ 
tion, with the meaning that the interaction is indepen¬ 
dent of the distance between the anyons and arises only 
from the inherent anyonic statistics. In Refs. 23 and 24 
the authors report on some phases that appear in lattice 
models of itinerant anyons, where the anyons—coupled 
by a Heisenberg interaction—are located on the sites of 
a lattice, with vacancies which allow for anyons to hop 
between sites but without braiding around one another. 
In Ref. 25 the authors study the real time dynamics of a 
single anyon moving between the sites of a ladder lattice 
with static anyons pinned to the plaquettes of the ladder, 


which serves as a model of coherent noise in topological 
quantum memories, and uncovers a signature that distin¬ 
guishes abelian anyons from non-Abelian anyons based 
on their transport properties. Noise models for medium 
sized topological memories based on real time stochastic 
dynamics of braiding Ising models anyons,^® Fibonacci 
anyons,and quantum double model anyons^® have also 
been studied. In this paper, we describe how to simu¬ 
late ground states of ID and quasi-lD models of itin¬ 
erant anyons, which may or may not involve braiding, 
and possibly include a Heisenberg interaction. We bench¬ 
mark our method by reporting ground state energies and 
ground state entanglement for these models. 

Large anyonic systems, like generic quantum many- 
body systems, are hard to simulate on a classical com¬ 
puter due to the exponential growth in the dimension of 
the state space with the number of particles. Until re¬ 
cently, numerical studies of anyons have primarily used 
exact diagonalization,^^’^®’^®^®^ which limits analysis to 
small system sizes and relies on finite-size scaling to ex¬ 
tract properties in the thermodynamic limit. A more 
successful approach uses tensor networks (TNs) which 
describes quantum many body states using a network 
of low rank tensors which can be contracted together to 
compute relevant quantities such as ground state energy, 
correlations, subsystem entropy, etc. One of the simplest 
tensor networks is the matrix product state (MPS) which 
forms the basis of highly successful algorithms, namely, 
the density matrix renormalization group (DMRG)®®^®® 
and the time-evolving block decimation (TEBD),®®^®® to 
simulate the ground state and dynamics of ID and quasi- 
ID quantum many-body systems. Exploiting translation 
invariance in TN states has allowed the study of systems 
directly at thermodynamic limit, circumventing the lim¬ 
itation on size encountered in exact diagonalization.®® 

Owing to their success for spin systems, tensor network 
algorithms have recently been adapted to simulate quan¬ 
tum many-body systems of anyons.In particular, 
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anyonic versions of the Matrix Product States (MPS), 
and of the TEBD and DMRG algorithms have been pro¬ 
posed and tested with a high degree of accuracy for any¬ 
onic chains.Tensor network algorithms are adapted 
to anyons by explicitly hardwiring the constraints im¬ 
plied by the fusion rules of the anyon model into the ten¬ 
sor network ansatz. This provides two important advan¬ 
tages. First, an anyonic TN representing a many-body 
anyonic state contains fewer complex coefficients than 
a non-symmetric TN description of the same state that 
does not explicitly encode the anyonic symmetry, thus 
providing for computational speedup. Secondly, using an 
anyonic TN as an ansatz in numerical simulations guar¬ 
antees that one remains in the physically relevant sector 
of the Hilbert state, namely, one with the desired to¬ 
tal anyonic charge, and thus avoiding leakage into states 
that are not allowed by the physics of the system, due to 
numerical errors. 

In this paper, we describe how to simulate the ground 
state of a system of itinerant anyons by means of the 
anyonic TEBD algorithm that additionally incorporates 
a U(l) symmetry corresponding to conservation of par¬ 
ticle number density. Our construction of the combined 
Anyon x U(l) symmetric MPS is the first to allow for 
simulating these systems with an arbitrary, specified ra¬ 
tional particle number density (or filling fraction), and 
gives direct access to Hilbert space sectors enumerated 
by anyonic charge and particle number density. Our 
MPS ansatz also allows us to simulate bosons, fermions, 
and anyons using the same algorithm, since bosons and 
fermions can be treated as simple types of anyons. 

Models of itinerant hardcore particles (bosons, 
fermions or anyons) in one dimension all have the same 
ground state properties since the particles do not ex¬ 
change positions. However, in two or higher dimensions, 
there are several paths by which particles may exchange 
positions. Therefore, beyond ID, ground state prop¬ 
erties of hardcore bosons, spinless fermions and hard¬ 
core anyons should reflect the influence of their exchange 
statistics. We test our method using itinerant Fibonacci 
anyons on a chain and itinerant braiding (henceforth, 
simply “braiding”) Fibonacci anyons on a two-leg lad¬ 
der, and show how the ground state energies differ from 
those of hardcore bosons and spinless fermions. We also 
present results for the ground state entanglement of the 
Golden Ladder model comprised of Fibonacci anyons in¬ 
teracting by means of ferromagnetic or antiferromagnetic 
Heisenberg interactions. 

Thematically, the paper is divided into two parts. The 
first part develops the Anyon x U(l) symmetric TN for¬ 
malism, and the second part describes applications of this 
ansatz to the simulation of models of itinerant and braid¬ 
ing anyons. The Anyon x U(l) symmetric MPS com¬ 
bines the recently proposed anyonic MPS'^^’"*’^ with the 
implementation of a U(l) symmetry in the MPS"*^®’^® and 
as such our presentation contains some review of both 
elements separately, which serves both as a reminder of 
important concepts and also introduces useful terminolo- 




FIG. 1. (i) Splitting tree and (ii) fusion tree, defining the “kef 
and “bra” bases respectively for a total number k of anyonic 
charges (ai, 02 ,..., Ufc) on the leaves, and (&i, 62 , •. •, hu- 2 ) as 
the fusion products on the links of the trees. The charge c 
on the trunk can, in principle, take all possible charge values 
permissible by the anyon model. If however, the fusion tree 
defines the basis of a pure quantum state, the charge c can 
only be the vacuum charge I. 


gies that persist throughout the paper. The structure of 
this paper is as follows: in Sec. H, we review the any¬ 
onic MPS. In Sec. HI, we review the implementation of 
a U(l) symmetry in the MPS corresponding to conser¬ 
vation of particle number density, in particular showing 
how it can be achieved as an instance of the anyonic MPS 
and how an arbitrary filling fraction is realized at the 
level of the ansatz. In Sec. IV we construct the combined 
MPS ansatz that incorporates both the anyonic symme¬ 
try and the U(I) symmetry. We present test models and 
benchmarking results in Sec. V and some conclusions in 
Sec. VI. 


II. ANYONIC MATRIX PRODUCT STATES 

We give a brief review of the anyonic MPS constructed 
in Refs. 43 and 44. For more details, the reader can 
consult these articles. 

The basis of the Hilbert space of anyons is described by 
a labeled directed fusion tree (see Fig. 1) where the charge 
c on any incoming edge at a vertex is determined from 
the charges a and h of the two outgoing edges around the 
same vertex, according to the fusion rules of the anyon 
theory 

a X b —>■ ^ c, ( 1 ) 

C 

which implies that charges a and b are allowed to fuse 
to, possibly, several different charges c. The is the 
multiplicity tensor, which encodes the number of ways 
of obtaining charge c from charges a and b. We con¬ 
sider only multiplicity-free anyon models in this work, 
with ATj, = 0 , 1 , which includes some of the models most 
relevant to current experiment such as Ising anyons and 
Fibonacci anyons. When ^ anyon model 

is non-Abelian. Anyonic charges have quantum dimen¬ 
sions analogous to the dimension of an irrep for a group, 
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and the dimensions da^ db-, and dc of three charges a, 6, 
and c must satisfy 

d,db = Y,N:bdc, ( 2 ) 

c 


absence of a charge. 

Formally, for a lattice £ of L sites with anyonic charges 
oil = (ai,Mai)jQ ^2 = ( 02 ,/iaa), ■ ■ ■ , = (aL,Mai,)i f^e 

anyonic MPS encoding the ground state 'Fgs is given 
diagrammatically as 


though in contrast with group theory, there is no require¬ 
ment that the quantum dimensions be integer. The total 
quantum dimension is then defined as I? = da sum¬ 

ming over all anyon charges a of the theory. 

The labeled fusion/splitting tree in Fig.l contains 
many charge labels, and can be extremely verbose when 
dealing with large anyonic systems. While explicit label¬ 
ing of fusion trees is possible in principle, it is not very 
practical for anyonic tensor network simulations. A bet¬ 
ter alternative is to enumerate the labeled fusion trees 
having a particular charge c at the trunk of the tree. To 
this end, let c is the total charge at the trunk of the fusion 
tree and introduce a new index Hc that enumerates each 
unique labeled fusion tree in increasing numerical order, 
/ic = 1, 2, • • • , z/c. Here is called the degeneracy index, 
and Vc is the degeneracy of the charge sector All the 
fusion trees are therefore concisely labeled by the multi¬ 
index 7 = (c, /ic)i with c as the total charge label and 
as its degeneracy index. A tensor network consists of con¬ 
nected tensors, which may be a combination of, single¬ 
index tensors (or vectors), two-index tensors (or matri¬ 
ces), or multi-index tensors. Similarly, the tensor objects 
of our Anyonic-U(l) MPS ansatz are the anyonic analog 
of the non-symmetric tensors, although, for anyons, we 
do not permit using tensors with more than three legs, 
as the fusion tree labelling of non-Abelian tensors can 
no longer be uniquely specified purely by multi-indices 
on the external legs. The details of how to construct 
anyonic tensors are given in the Appendix A. 

One convenient form of the conventional MPS ansatz 
is that given by Vidal,which is an array of two-index 
and three-index tensors forming a linear network of ten¬ 
sors. For a finite lattice with open boundary condition, 
the tensors on the boundary of the MPS (i.e. the first and 
last sites) are two-index tensors while the “bulk” of the 
network consists both of two-index tensors (Schmidt vec¬ 
tors) and three-index tensors for each of the other {n — 2) 
sites. 

Analogously, the MPS was adapted to anyons by Singh 
et. al. in Ref. 43, using the basic anyonic tensors (two- 
index and three-index anyonic tensors) after the pattern 
of the conventional MPS. Each three-index tensor is in¬ 
dexed by both the charge and the degeneracy of the 
anyons making up each site. The charges on the trivalent 
vertex of the tensor are compatible in accordance with 
the fusion rules of the anyon model. The Schmidt vectors, 
which are two-index tensors, are charge-conserving diag¬ 
onal matrices. The basis labeling ai = (oij/iaJ for each 
site of the anyonic lattice is given by the set of charges 
ai and the degeneracies /Xq. of each charge. The labels 
piai take fixed value 1 if there is only one configuration 
for each possible charge labeling at each site, e.g. if the 
possible physical states are merely the presence or the 
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are obtained by an iterative fusion of the multi-indices 

[di ipi^ fibi) ^-od Oj_|_3 (^i-t-17 Mai+i )7 

Pi X Oi+i —>• /3i+l7 (4) 


where, as before, the charge 6i+i, 

h+i= Y. <:i,(&*xa,+i), (5) 


and the total degeneracy ^bi+i of the charge 6i+i is de¬ 
termined by 


Vb,+i = 


bi,ai+i 


K 




( 6 ) 


It should be noted that this anyonic MPS has been 
drawn with site indices going upwards, to make apparent 
the visual similarity with anyonic fusion tree diagrams, 
but it is essentially the same ansatz given in Ref. 43. Due 
to the iterative fusion process down the tree of the any¬ 
onic MPS the dimensions of the tensors Pl'l required to 
exactly construct an arbitrary state will vary, but in prac¬ 
tice an upper bound is imposed on the bond dimension 
X ahead of time. The bound chosen usually depends on 
the amount of entanglement and correlations needed to 
faithfully represent the state of the system (and on com¬ 
putational resources available). As such, anyonic MPS 
provides a systematic way of handling anyonic systems, 
specifying both the basis (i.e. the fusion tree) and encod¬ 
ing the amplitudes of the state in the tensors. 

As a proof-of-principle example, this anyonic-MPS 
ansatz has been used to simulate, together with the 
anyonic-TEBD algorithm, a chain of interacting non- 
Abelian anyons (e.g. Fibonacci and Ising anyons) cou¬ 
pled by a Heisenberg interaction. The charge multi-index 
ai on each site i of the leaves of the anyonic-MPS is set 
(in the case of Fibonacci anyons) to ai = (r, 1), where t 
is the Fibonacci anyon charge, and the number 1 is the 
degeneracy of the r charge on site i (i.e. the number of 
different configurations on the site consistent with a to¬ 
tal charge of r). The anyonic MPS is, however, a general 
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ansatz capable of dealing with systems with any quantum 
group symmetry, and hence, can be adapted to work with 
other symmetries. Abelian or non-Abelian. For instance, 
by replacing the anyonic charges with particle number 
charges, the anyonic-MPS can serve as a U(1)-MPS^®, 
which can be used to simulate physical systems having a 
global particle number TV on a finite lattice L. 

On an infinite lattice with translation invariance of the 
Hamiltonian, if the U(l) charge is identified with particle 
number then the U(1)-MPS is primitively a zero-density 
ansatz [i.e. one favouring a mean U(l) charge per site 
of 0], and cannot directly be used to simulate an infinite 
lattice with a finite non-zero particle density. In the next 
section we show how to tune the U(1)-MPS to simulate 
an inhnite lattice system at non-zero density, and in Sec¬ 
tion IV we propose a modified ansatz, the Anyonic-U(I) 
MPS, that conserves both particle density and anyonic 
charge symmetry, and which can be used to simulate any¬ 
onic systems (including braiding of anyons) at a specified 
rational filling fraction. 

III. U(1)-MPS AND PARTICLE DENSITY 
CONSERVATION 

In the last section we alluded to the fact that the any¬ 
onic MPS can serve as a U(1)-MPS by replacing the any¬ 
onic charge labels with the particle number charge labels. 
Specifically, let us consider a lattice £ of L sites, where 
each site can accommodate a finite number of particles, 
n = 0,1, 2,..., d — 1. The positive integers n can be re¬ 
garded as the irreps of the 17(1) symmetry, which can 
intuitively be understood as: n = 0 is the absence of a 
particle, n = 1 is the presence of one particle, n = 2 is the 
presence of two particles, and so on. The total number 
of particles N on the lattice of L sites is V = 
with a particle density of v = N/L. 

The Hilbert space of the lattice, can 

be alternative written as, ^ direct sum 

over subspaces with fixed numbers of particles n. Utiliz¬ 
ing this alternative structure a particle-number conserv¬ 
ing Hamiltonian H can be directly diagonalized in the 
Yn subspace, offering savings on the computational cost. 
The U(1)-MPS ansatz for N particles on an L-site lat¬ 
tice can be derived from the anyonic MPS by fixing the 
particle number N and degeneracy = 1 at the “right 
end” of the last tensor, and charge 0 (i.e. zero) on the 
“left end” of the first tensor. The on-site multi-indices 
of the “bulk” [L — 2) tensors carry where 

Hi is the U(l) charge on site i, and fii enumerates the 
degeneracy of that charge, for alH G £. The MPS bonds 
also carry charge and degeneracy indices, but unlike sys¬ 
tems of anyons where degeneracy comes from the fusion 
rules of the anyon model, degeneracy in U(l)-symmetric 
lattice models comes from the number of combinatorial 
arrangement of the charges on the lattice. 

Therefore, with a properly constructed ansatz and an 
optimization algorithm like TEBD or DMRG, one can 
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FIG. 2. An example of an infinite MPS with a block made up 
of three tensors F^^^ and F^^b By translational invariance 
of the Hamiltonian, an infinite MPS corresponds to an infinite 
repetition of the block and hence optimization to ground state 
is performed only on the tensors within a single block. 

compute the ground state of a local U(l)-symmetric 
Hamiltonian on a finite lattice. If this finite U(I)-MPS 
is naively extended to simulate an infinite lattice model, 
the ansatz would correspond to a zero-density ansatz be¬ 
cause of the finite size of the bond dimension y and the 
assumption that the U(l) charge labels exhibited on this 
bond are finite. In the next subsections we give a heuris¬ 
tic proof of this statement, and we then propose a tech¬ 
nique which can be employed to tune the U(1)-MPS away 
from being a zero-density ansatz, to any desired non-zero 
particle density. 

A. Zero-density U(1)-MPS 

Restricted to a hnite bond dimension y carrying finite 
U(l) charges, the U(1)-MPS with integer charge labels on 
an infinite lattice is a zero-density MPS ansatz. Consider 
a section of the infinite MPS in Fig. 2 with the charge- 
degeneracy indices on physical sites 

— (ui ,), 02 — (^ 25 /^ 712 ) 7^3 — (^ 37 /^ 713)7 

and on the links 

/3l = (mi,/XmJ,/32 = {m2,fim2),P3 = 

Pi = {mi,fj,,ni)- 

The on-site charges are set to take positive integer 
charges corresponding to particle number (e.g. hardcore 
boson has ni€ {0,1}). The charges rrii on the links take 
only a finite number of charges with degeneracy index 
/imi = 172 ,-•• ,Vmi- The charges and degeneracies on 
the bond are constrained by the finite bond dimension y 
and given as y = X)m where labels the charge 
on the link i. For any realistic computer simulation, the 
charge labels on the MPS bonds are all finite. Assume we 
cut the infinite lattice into two partitions. There exists a 
finite amount of charge k on the link of the left partition, 
corresponding to a finite number of particles, and the 
density on the left half-chain is therefore z/ = fe /00 —>■ 
0 and therefore the infinite U(l) MPS is a zero-density 
ansatz. However it is possible to remedy this and have a 
nonzero density U(l) MPS by shifting the on-site charges 
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so that a U(l) charge of zero corresponds to the desired 
filling fraction. We present this transformation below. 


B. Non-zero density U(1)-MPS 


By employing translation invariance, an infinite U(l)- 
symmetric MPS consists of a block of repeated U(l)- 
symmetric tensors, albeit that such an ansatz is zero- 
density and will yield a ground state of an empty lat¬ 
tice as seen above. However, by transforming the on-site 
charges of the MPS, we can cause a U(l) charge of zero 
to correspond to the desired density. 

For simplicity and without loss of generality we con¬ 
sider hardcore particles, with charge labels nG {0,1} on 
each site of the U(1)-MPS lattice. Let the desired den¬ 
sity on the infinite lattice he v = p/q, which can be 
interpreted as having an average of p particles on ev¬ 
ery q sites. Using the additive (abelian) fusion rules of 
U(l) charges, a U(1)-MPS with p particles corresponds 
to having p sites with charge n = 1 and the remaining 
q — p sites with holes n = 0. In an infinitely increasing 
block, the number of particles p increases infinitely, but 
by “subtracting off” the p number of particles, we can re- 
center the relevant subspace to be labeled by the charge 
0, which is retained in a practical simulation. Formally, 
by using the transformation. 


/ 

n 



= qn-p, 


P <q, 


the on-site charges transform as, 


n = 0 n' = —p 

n = 1 n = q — p 


( 7 ) 


where multiplication by q in Eq. (7) is purely for conve¬ 
nience and ensures that the n' charges, like the n charges, 
are integer. In essence, before this transform, the desired 
filling fraction in the MPS would correspond to having 
p occurrences of charge 1 and q — p holes 0, summing to 
a total charge of p. But after the transform, the desired 
filling corresponds to having p occurrences of particles 
with charge q — p, and q — p holes with charge —p, which 
sums to a total charge of zero. The charge distribution 
on any link on a U(l)-symmetric infinite MPS is cen¬ 
tered on the zero charge sector, which now corresponds 
to a particle density of p/q. Thus it becomes possible to 
tune the U(1)-MPS to the desired filling fraction without 
using tensors with more than three legs.^® An example of 
how this transform applies to the half-filling is presented 
as an example below. 


1. Example: Half-filled MPS ansatz for hardcore bosons 

Consider a particular configuration of an infinite lattice 
at half filling, where there is on average, one particle on 
every two sites as shown in Fig. 3(a). Each box represents 


0 

+1 

+1 
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+1 
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( Original charge labels ) 




-1 

+1 

+1 
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+1 
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( “Shifted” charge labels ) 


FIG. 3. Schematic representation of an infinite lattice with a) 
a typical half-filled configuration with one particle on every 
two sites, and b) a “shifted” version of a) but with average of 
zero particle on every two sites. 


a site and the charge on the site is indicated inside the 
box. There is on average one particle for every two sites, 
and assuming that this average density is maintained, 
this will correspond to half filling on the infinite lattice. 
This is of course not the only way to achieve half filling, 
but the example will suffice to illustrate how to achieve 
a half-filled U(1)-MPS. 

With only nonnegative charges on each site, i.e. 
nG {0,1}, the charges on the links of the MPS—which 
are derived by fusion of all charges leading to that link— 
are also all nonnegative. However the implementation 
requirement that the charge indices be finite (and the fi¬ 
nite size of the bond dimension) places an upper bound 
on the set of charges on the links which are retained af¬ 
ter truncation of the Hilbert space of the link. Hence the 
dominant larger-iV states in the infinite lattice are trun¬ 
cated. However, by using n' = 2n — 1, the on-site charges 
are re-defined as 0 ^ — 1 and 1 —s- 1, to give the “shifted” 
configuration in Fig. 3(b), for which the dominant states 
now inhabit the zero particle sector. Nearby charge sec¬ 
tors such as ±1 on the bonds represent small fluctuations 
in filling fraction relative to a baseline oi v = 1/2. We 
emphasize that the complex amplitudes of the state are 
not changed, only that their index is relabeled. 


IV. ANYON X U(1)-SYMMETRIC MPS 
A. Composite charges and fusion rules 

In the last Section, we reviewed the U(l) MPS and ex¬ 
plained how to achieve an arbitrary rational filling frac¬ 
tion on the infinite lattice. In this Section, we investigate 
how anyonic systems at arbitrary filling fractions can be 
simulated using an ansatz that conserves both the any¬ 
onic (quantum group) symmetry and the U(l) symmetry. 

We first recognize that the two symmetry groups 
are described by particle spectra with differing fusion 
rules. Similar to creating a new group from product 
of two groups, we introduce the Cartesian product of 
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FIG. 4. Fusion of composite charges situated on a manifold 
which supports either a single nontrivial anyonic Fibonacci 
charge (r, 1) or a vacuum charges (1,0) at each point. A 
lattice may then be embedded into the manifold. In the 
figure, a linear ordering has been defined and fusion pro¬ 
ceeds from left to right. The total charge outcomes which 
are ((I, 3), (r, 3), (t, 3)), indicates that there are in total three 
Fibonacci charge r on the manifold fusing either into the vac¬ 
uum charge (I) or the Fibonacci anyon r channel. As there 
are two charges with label (t, 3), we would also introduce a 
degeneracy index fJ,(T, 3 ) = 1, 2 to enumerate these outcomes. 


along with the composite fusion rules. To have an ansatz 
for a particular anyonic filling fraction, the method of 
shifting the U(l) charges can be employed. This only 
amounts to a shift in the U(l) charge labels, while the 
labels on the anyonic fusion space are not altered. The di¬ 
agrammatic representation of tensors with the new sym¬ 
metry group and the MPS ansatz constructed from them 
are the same as given in Section A and we will not repro¬ 
duce them here. 

B. Manipulations of Anyon x U(l) tensors 

Topological manipulations such as F-moves, R-moves, 
and vertical bends applied to anyonic fusion trees are also 
modihed in the case of a Anyon x U(l) symmetry. Let 
the label a = (a, n) be the composite charge where a is 
the anyonic charge and n is the U(l) charge. Below we 
present typical manipulations needed to contract anyonic 
tensors during optimization of anyonic MPS. 


the anyonic charge spectrum A = {a, ,d} and 

the U(l) charge spectrum which will be designated as 
U = {n, m, • • • , z} where n are integer charges, n G Zoo. 
The product of the two particle spectra is given as 
AxU = {(a, n) \ a G A,n G U}, where the label (a, n) 
is referred to as the composite charge. The charges on 
the physical site and on the links of the MPS are taken 
from this set AxU. 

The “new” fusion rules for the composite charges are 
derived from the fusion rules of the two theories. 


(ai,ni) X (0,712) = (oi X 02,711 x 712) 

= (^ma 2 «i 2 , ni + 712) , ( 8 ) 

012 

where as aforementioned rii x has a unique outcome 
{rii + 712 ) with an additive fusion rule, while the non- 
abelian anyons have generally more than one fusion out¬ 
come, hence the need for the summation Y” over all 
possible charge outcomes 012 . 

We consider only hardcore anyonic particles, meaning 
that either there is a nontrivial anyonic particle on a site 
or the site is vacant. The vacuum charge of the composite 
charge spectrum AxU is (I, 0). The presence of a single 
nontrivial anyonic charge is represented by (a, 1) where 
a e M \ I and the U(l) charge 1 imposes a hardcore 
constraint of a single charge on the site. The use of the 
U(l) charge allows the counting of the anyonic charges 
fusing into a particular fusion channel irrespective of the 
outcome anyonic charge. A simple example is shown in 
Fig. 4. 

The anyonic MPS ansatz and the U(l) symmetry dis¬ 
cussed in previous Sections can be used together to re¬ 
alize an Anyon x U(l)-syinmetric MPS ansatz with the 
desired particle density. The minor modification needed 
in the new ansatz involves using the composite charges 


1. F-moves 


The first topological manipulation required is that of 
changing the fusion order of the composite charges repre¬ 
sented by the fusion tree. Let the basis fusion tree where 
fusion of charges proceeds from left to right be referred 
to as the standard basis. If instead a different fusion or¬ 
dering is chosen, such as fusion from right to left, the 
charge outcomes are still the same, a fact guaranteed by 
the constraint of associativity. Formally, this associativ¬ 
ity constraint corresponds to the Pentagon Equations, 
as given in e.g. Ref. 49. The corresponding operation if 
the F-move, which transforms from one fusion basis to 
another one and is given diagraminatically as 

d b c d b c 

=E/(7|“)' 

d d 

(9) 

where the coefficient decomposes into its any¬ 

onic and U(l) counterparts as, 

{Fr)[ = {Fsni (w) 


The factor is given by the F coefficients of 

the anyon model while the U(l) factor is given by 
^ which equals 

one if the charges are compatible or zero otherwise. F- 
moves may also be applied to pairs of contiguous vertices 
appearing within a larger diagram. 

It was noted in Ref. 50 that a symmetric tensor de¬ 
composes into a linear superposition of the degeneracy 
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tensor and its spin network for systems with nontrivial 
symmetries such as SU(2), and more generally also for 
quantum symmetries. Therefore any section of the any- 
onic MPS can be decomposed into its degeneracy tensor 
and anyonic network as 



a 


b 

fd'b 

C 

/^c 

0 

1 

0 

1 

0 

1 

0 

1 

1 

1 

1 

1 

1 

1 

0 

1 

1 

2 

1 

1 

1 

1 

0 

2 


TABLE I. The mapping from tensor product state \a,^a) (8> 
\b,iJ,b) to a new basis \c, ^c) using the Z 2 fusion rule. Degen¬ 
eracy basis labels for each charge x £ (a, b, c) have been 
included to count fusion into a particular charge. 


The F-move is then applied on the anyonic diagram and 
the resulting F-factors are absorbed into the tensor re¬ 
sulting from contraction of the degeneracy tensor net¬ 
work. As shown, this process is valid for any portion of 
the diagram where the F-move operation can be applied. 


2. R-moves 

Anyons have very rich particle exchange statistics 
which are neither bosonic nor fermionic. The exchange 
factors are encoded in the R-matrix which is a matrix rep¬ 
resentation of the braid (or R-) move. The braid operator 
for composite anyonic charges is given diagrammatically 
as 



dim(V*^‘^^) = dim(V(^)) x dim(V(^)). The linear map 
can be written as 

T=J2 Ta, \c) {a\ 0 {b\, (12) 

a,6,c 

which sends a product basis |a) 0 |6) to the basis state 
|c). This assignation can be manually defined by, for in¬ 
stance, iterating slowly over the basis labeled by a and 
fast over the basis labeled by b, sending them to a unique 
new basis indexed by c. The coefficients of are 1 for 
a valid map {ab —> c) and zero when there is no valid 
map. Consider the following example. Let be a 

vector space of dimension = spanc{|a;)^}^lQ ^ 

and similarly let be a vector space of dimension 

ds' = spanc{|?/)g}(j®Q The tensor product 

space 0 can be labeled by states 

{|z)^} using the simple assignment map |x)^ 0 \y)g —> 
\z = dsx + y)(.. 

However, more structure can be included into the lin¬ 
ear map by defining some relationship between the basis 
labels of the spaces. For example, assume we include a 
Z 2 fusion rule defined by 

1x1—>-0; Oxa^-a Vo. (13) 


where the factor decomposes as 

_ j^ab 

and = 1 if no + n& = ric. The factors are given 

by the anyon model. 

To model the braiding of anyons by exchanging the 
positions of anyons, the Hamiltonian H should contain 
the braid operator. Later, we construct a Hamiltonian 
for the braiding of anyons supported on the vertices of a 
ladder. 


3. Fusion tensor and loop factors 

A trivalent tensor can be used to define a linear map 
from the tensor product of two Hilbert spaces and 
(which can possibly be degenerate) to a new com¬ 
posite space The dimension of the new space 


The charge outcome c = 1, resulting from the fusion 0x1 
and 1 x 0, is degenerate, as is c = 0 which results from 
0x0 and 1x1. The degenerate outcomes are then indexed 
by a degeneracy index The linear map using the Z 2 
fusion rule is given in Table 1. 

Therefore, as per Ref. 45, the linear map tensor can in 
general be written as 

E iT^btU\c,pi,){a,ya\{b,yb\, (14) 

a,b,c 

where the tensor T is constructed blockwise from ten¬ 
sors with each block being identified by the charge 
triple {a,b,c). Each block tensor then has its en¬ 
tries indexed by the corresponding degeneracy indices 
(/^aj M&I Re)- 

We generalize this to anyonic systems admitting 
Anyon x U(l) symmetries as follows: Let two sites of 
an anyonic system be described by a degenerate Hilbert 
space and with basis {a = (d,/ra)} and 








a) 



FIG. 5. a) The diagrammatic representation of the anyonic 
fusion tensor N which can be expressed in its block structure 

^ (^) ^ map a X ^ 7 and 

zero otherwise . b) The fusion tensor and its Hermitian 
conjugate iV“,^ are linear and hence their product contracts 
to the identity operator defined on the new multi-index 7 . 
c) Elimination of loops from anyonic diagrams as given in 
Ref. 49. 


{/? = {b, and let the anyonic fusion product define a 
“fusion map” iV^ ^ from multi-indices a and /3 to a new 
multi-index 7 . The anyonic fusion map creates a new ver¬ 
tex and we normalize it according to diagrammatic iso- 
topy convention. As was discussed previously, the map 
is created by iterating slowly over basis label a and fast 
over /3, and enumerating pairs (a, /3) by a new label 7 . 
The fusion tensor is represented in Fig. 5(a). However, 
unlike the case of abelian symmetry, for anyons normal¬ 
ized according to the diagrammatic isotopic convention 
the coefficients of a valid fusion map a x /3 —>■ 7 take the 

/ , \i/4 

value of the vertex normalization factor ■ As 

for abelian anyons, the coefficients are zero if there is no 
valid fusion map. 

The anyonic fusion tensor and its Hermitian con¬ 
jugate, the splitting tensor are linear maps and 

fulfill the condition that N^p = ly (Einstein sum¬ 
mation convention assumed) which is an identity oper¬ 
ator on the new (degenerate) space as shown in 

Fig. 5(b). The loop resulting from the contraction in 
Fig. 5(b) is eliminated using the relation in Fig. 5(c). It 

should be noted that the vertex normalization in 

the definition of the fusion tensor N and splitting ten¬ 
sor iVl cancels with the loop factors and hence 

the identity matrix operator in Fig. 5(b) does not con¬ 
tain any factor of the quantum dimension da of anyonic 


charge d. Also note that the quantum dimension da de¬ 
composes as the product da = dadn^ where da is the 
anyon quantum dimension and dn^ is the dimension of 
U(l) charge, which is trivially equal to one. 

4- Vertical Bends 

Bending a charge line horizontally is trivial, as time¬ 
like (i.e. horizontal) slices of the fusion tree are invariant 
under topology-preserving deformations. However, ver¬ 
tically bending an anyon charge line is non-trivial and 
involves reversing the orientation of the anyon worldline. 
The details of how to resolve the vertical bends in terms 
of F-moves have been given in Ref. 49 and also in Ref. 44. 
We do not repeat the derivations here but only mention 
the minor changes in the presence of U(l) charges. 

We adapt the left bend given in Ref. 44 to the case of 
Anyon X U(l). This is given as 



(15) 

where the U(l) charges on the vertex satisfy the condition 
na + rij, = ric- The dual of the anyonic a and that of 
the U(l) charge n are denoted respectively as d and — n, 
which will be the new charge label on the left-bent leg. 
In addition there is also an implicit F coefficient from 
the U(l) charge sector, but this is always equal to 1. 
Similarly, the right bend is given by 



In summary, by constructing the appropriate tensor 
objects with Anyon x U(l) symmetry (e.g. two index 
and three index tensors, fusion tensors, etc.), one can 
construct an MPS ansatz for anyonic systems at any 
fixed rational filling. This ansatz may then be used to 
construct an an approximation to the ground state of a 
system by means of anyonic algorithms such as the any¬ 
onic TEBD algorithm proposed by Singh et al.^^ or the 
anyonic DMRG.'^^ Ground state properties such as en¬ 
tropy scaling and correlation functions can be computed 
using approaches similar to those for conventional ten¬ 
sor networks, but modified to account for anyonic statis¬ 
tics by normalizing vertices, removing loops and bending 
anyonic charge lines in accordance with the prescriptions 
given in Ref. 49. 


V. TEST MODELS AND RESULTS 

We provide some examples to demonstrate that an 
Anyon x U(l)-symmetric MPS ansatz may be used to 
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simulate itinerant anyonic systems at any rational fill¬ 
ing fraction, and also to provide an example of a ten¬ 
sor network where more than one symmetry is exploited 
in the algorithm, with one of these symmetries being 
anyonic. The anyonic models we consider are itinerant 
hardcore Fibonacci anyons with variable filling fractions 
and interactions, on a one-dimensional chain (the Golden 
Chain^®), and on a ladder (the Golden Ladder). We 
compute their ground state energies and entanglement 
entropies, using the definition of entanglement entropy 
for non-Abelian anyons given in Ref. 51. From this we 
extract the central charges of the conformal field theo¬ 
ries associated with the infra-red limits of these models. 
Analytical solutions for these models are not generally 
known, but we establish the validity of our method by 
using it to compute equivalent known results for spinless 
fermions and hardcore bosons, and also by comparing 
results for selected anyonic systems with those obtained 
using anyonic DMRG."’^'* In general, our results are found 
to be accurate to 4 or 5 decimal places. 


means has two eigenvalues, 0 (when s = 0) and 2 
(when s = 1). Therefore can be written in terms 
of projectors to the singlet and triplet subspaces as 
(Si + Sj)'^ = where and are the 

projectors to singlet and triplet subspaces. Therefore, 

4 + i, (21) 

where the identity, I = -|- has been used in the 

last step. Therefore, the t-J Hamiltonian simplifies to 

H = —f ^ J ^ + const. (22) 

<b> <b> 

For J > 0 the Hamiltonian favours neigbouring spins 
forming singlets (antiferromagnetic), and for J < 0 it 
favours triplet formation (ferromagnetic). We adapt the 
electronic t-J to anyons. 


A. Itinerant hardcore particles on a 
one-dimensional chain 

We give some diagnostic test results for hopping and 
interacting anyons on a chain using an anyonic t-J Hamil¬ 
tonian which is analogous to the t-J model for electrons. 
To make the analogy more apparent, we briefly review 
the electronic t-J model. 


1. Electronic t-J Model 


The electronic t-J Hamiltonian consists of two com¬ 
peting terms: a term corresponding to the kinetic energy 
of the electrons, and an interaction between their spin 
degrees of freedom. The t-J Hamiltonian is 

H = -t^c\cjS jY^Si-Sj, (17) 

(b> (b> 


where the first term is the kinetic energy with hopping 
strength t and c|(ci) is the creation (annihilation) oper¬ 
ator which satisfies fermionic anticoinmutation relations. 
The second term is the Heisenberg spin-spin interaction 
which can be rewritten in terms of projector of nearest 
spins to the singlet state using the fact that 


■Sj = - 
2 


(S. + 


(18) 


The addition of two spin-1/2 charges is given by the rule. 


1 

2 


1 

G - = 00 1. 


(19) 


Let S = Si-\- Sj and choose units such that h= 1. Then 
the relation 


5^ |s, m) = s(s + 1) |s, m), (20) 


2. Anyonic t-J Hamiltonian in ID: Hopping term 

Anyonic operators are written as matrices on the fu¬ 
sion space of the participating anyons. A local two-site 
Hamiltonian H = ■*■^1 can be written diagram- 

matically as 





(23) 

where the values of the function i7(ai, a^+i, o', c) 
are determined by the model being constructed. The 
conservation of charge c resulting from fusion puts the 
Hamiltonian into block-diagonal form as = 

The anyonic t-J Hamitonian consists of two terms: a 
hopping term and an interaction term. To give a system¬ 
atic and concrete treatment of both terms. We give the 
explicit construction for Fibonacci anyons. 

The hopping of a Fibonacci anyon in ID means the 
neigbouring site has to be vacant, corresponding to the 
vacuum charge I. The kinetic operator can thus be rep¬ 
resented as 


D T T D 



r D Dr 


(24) 

being analogous to the fermionic terms (c|^Ci+i_o. + h.c) 
that translate a fermion between sites i and z + 1. Since 
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FIG. 6. (Color online) Ground state energy (E) of itinerant 
hard-core Fibonacci anyons on an infinite chain at different 
filling fractions {v). The data points resnlt from numerical 
simulations, while the smooth cnrve is the gronnd state en¬ 
ergy for an equivalent system of spinless fermions at the ther¬ 
modynamic limit. 

the anyonic hopping term requires there is a vacant site 
with vacuum charge I, this implies that the hopping term 
is nonzero only when (a^ = I, a^+i = r, a' = t, = 
I) or when (a^ = t, a^+i = I, a' = I, = r). For 
dynamics in one dimension with hardcore constraints, the 
underlying exchange statistics of the particle do not affect 
the ground state properties (though the degeneracy of the 
ground states may differ for different particle species). 
Therefore, itinerant hardcore Fibonacci anyons, spinless 
fermions and hardcore bosons all have the same ground 
state energies at any rational filling of the lattice. 

We use the Anyon x U(l) symmetric TEBD algorithm 
to compute the ground state energies of itinerant Fi¬ 
bonacci anyons, spinless fermions and hardcore bosons 
on a ID lattice. We obtained the same ground state en¬ 
ergies for these three cases up to 4 to 5 decimal places. 
This is owing to the fact that particles are not allowed 
to exchange positions on the lattice, and thus particle 
statistics do not affect the ground state propoeties. 

In Fig. 6 we plot the numerical ground state energy of 
itinerant Fibonacci anyons against the analytical ground 
state energy of an equivalent system of spinless fermions 
given by 

= (25) 

TT 

Numerical ground state energies for spinless fermions and 
hardcore bosons result in an identical curve and so are 
not reproduced here. 


3. Anyonic t-J Hamiltonian in ID : Heisenberg interaction 

term 

Next we include an anyonic Heisenberg interaction 
term in addition to the hopping term. The anyonic 


Heisenberg interaction is constructed by analogy to the 
Heisenberg spin-spin interaction. For 100% filling this 
model was first proposed and studied by Feiguin et al.,^® 
and is known as the Golden Chain. The anyonic Heisen¬ 
berg interaction takes the form 


T T T T 



T T T T 


where Jj > Jr corresponds to an antiferromagnetic in¬ 
teraction favouring fusion of the two Fibonacci anyons 
to the vacuum charge I, and Ji < Jr corresponds to a 
ferromagnetic interaction favouring projection to the Fi¬ 
bonacci charge r. 

When a Heisenberg interaction is introduced into a sys¬ 
tem of itinerant Fibonacci anyons, the extensive degener¬ 
acy of the free anyon system is lifted. The Hilbert space 
of the interacting itinerant anyon system admits the de¬ 
composition 

TJ — "^config ^ ^fusion (27) 

where F^config is the space of particle configurations, and 
^fusion is the space of valid labelings of the fusion tree. 
The Hamiltonian admits an equivalent decomposition, 
and the Hamiltonian for a system of free particles (act¬ 
ing on Hconfig is associated with a central charge of 1. 
When a Heisenberg-type interaction is added, this acts 
on Hfusion 7 lifting the degeneracy of the states in this sub¬ 
space. For a critical interaction, the total central charge 
is additive, and may be written 1-l-c where 1 is the contri¬ 
bution from the itinerant anyon model acting on F^conflg 
and c is the contribution from the interactions on the fu¬ 
sion portion of the Hilbert space.®® This is alluded to as 
spin charge separation. From our numerical simulations, 
when Jf > Jr (antiferromagnetic), we obtained c = 0.708 
and when Ji < Jr (ferromagnetic), we obtained c = 0.84, 
for total central charges of 1.708 and 1.84. These are very 
close to the expected central charges of 1 -1-7/10 for anti¬ 
ferromagnetic interaction and 1 -I- 4/5 for ferromagnetic 
interaction. 


B. Anyonic t-J Model on ladder 

Nonabelian anyons have nontrivial braid factors, mak¬ 
ing their simulation difficult. For such systems, numerical 
approaches based on Monte Carlo schemes are plagued 
by a form of the sign problem. We offer numerical evi¬ 
dence that anyonic tensor networks (such as the anyonic 
MPS) are able to simulate anyonic systems on geometries 
beyond one dimension, in situations where the anyons 
experience braiding. To model how braiding statistics 
affect the ground state of anyons, we introduce the any¬ 
onic t-J model on the ladder, as a generalisation of the 
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Jv 


FIG. 7. A ladder of itinerant anyons. Ficticious strings are 
attached to each nontrivial charge to indicate that they can 
participate in nontrivial braids as they exchange positions 
with neigbouring anyons. For example, the nontrivial any- 
onic charge Oi braids with the nontrivial anyonic charge &i_i 
as it hops horizontally to a new site with Oi-i = 0. The la¬ 
bels Jv and Jh represent the amplitudes for projecting the 
corresponding circled pairs into the vacuum sector. 



7 7 


FIG. 8 . (Color online) The two convenient fusion orderings, 
with their respective fusion trees shown underneath. The first 
fusion order couples charges (a, b) and (c, d) while the second 
couples charges (a,c) and {b,d). 


model already considered on a chain. Each site on the 
ladder supports only two types of charges, namely, either 
a vacuum charge I or a single Fibonacci anyon r. Un¬ 
like in one dimension, anyons on the ladder can exchange 
positions and consequently braid. 

To model the braiding of anyons on the ladder in a 
consistent manner, we impose a linear ordering to the 
anyons by attaching ficticious “strings” to the anyons 
and oriented them leftward of their on-site position (see 
Fig. 7). When an anyon hops from one site to another 
on either the top or bottom chain of the ladder, it braids 
with any adjacent anyonic charge along its trajectory, 
with the strings acting as a convenient mnemonic to vi¬ 
sualise the orientation of the braid. 

With reference to Fig. 7, the anyonic t-J Hamiltonian 
can be written as 

N-l 

H — tji ^ ^ I T “b h-C^ 

i=l 

- Jh (n" 

t ^ 

“ y X] + f>ai+i^bi+i=I + h-C^ 

^ i=l 

-Y(nL,..+nL+,.(28) 

where {th,tv) and {Jh,Jv) are the hopping and interac¬ 
tion amplitudes for anyons on the legs and rungs of the 
ladder. The vacuum charge is denoted by I. The oper¬ 
ator hx^y=i moves a nontrivial charge x into a new site 
having trivial vacuum charge y = I while it braids the 
charge x with any other charge along its path. The pro¬ 
jector n® y projects the nontrivial anyonic charges x and 
y into a vacuum charge I. The anyonic interaction is an¬ 
tiferromagnetic when J > 0 and the interaction becomes 
ferromagnetic when J < 0. The Hamiltonian along the 
rung has been symmetrized with half a contribution from 
each of the rungs on sites i and i -|- 1 . 


Below, we show an explicit derivation of the Hamilto¬ 
nian terms which can be arranged as a charge-conserving 
matrix operator. The local Hamiltonian h is derived on a 
plaquette whose vertices are labeled (a, 6 , c, d) for brevity 
as shown in Fig. 8 . The local Hamiltonian is written as 

h = —tfi -|- h-c^ — Jh T 

- I {ba^b=l + k^d=I + h-c) - ^ (fllh + . 

(29) 

Depending on the imposed fusion order, some of the 
operators will be diagonal in the fusion basis. The two 
most convenient fusion order are shown in Fig. 8 . Let 
the first basis be denoted as |I) = \{ab;a){cd; 
with the fusion order ((a, 6 )(c, d)) where the anyons 
(a, b) and (c, d) are first fused independently, then fuse 
their outcomes and let the second basis be |H) = 
|(ac; K)( 6 d; A)(kA; 7 )) with fusion order ((a, c)(b, d)). Us¬ 
ing a series of F-moves and i?-moves, the first basis trans¬ 
forms into the second basis according to 


abed a c b d 



7 7 (30) 


where the tensor is given by. 



(31) 


with its derivation being given in Appendix C 1. 

The rest of the derivation of the matrix expression for 
the Hamiltonian h is performed as an illustrative example 
in Appendix C2. 
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FIG. 9. (Color online) The ground state energies of hard¬ 
core bosons (HCB), spinless fermions and itinerant Fibonacci 
anyons (HCFib) on a two-leg ladder when only hopping is 
turned on. The line is a guide to the eye. 


1. Itinerant Fibonacci anyons, spinless fermions and 
hardcore bosons on a ladder 

We exploit the anyonic and U(l) symmetries of the 
model both in the MPS ansatz and in the Hamiltonian 
H, and use the TEBD algorithm to compute the ground 
state energies of itinerant Fibonacci anyons on the ladder 
at different filling fractions. Since the MPS has a one¬ 
dimensional structure, we map the ladder to a chain by 
fusing the anyonic charges on each rung to make a new 
single site. The vertical and the horizontal hopping rates 
are set equal to one, th = 1 and ty = 1 while the vertical 
and horizontal Heisenberg interactions Jy and Jh are set 
to zero. 

There are no known analytical results for the ground 
state of itinerant Fibonacci anyons on a ladder, but 
we test the validity of our method against the ground 
state energies of itinerant hardcore bosons and spinless 
fermions on the ladder shown in Fig. 9. The phase dia¬ 
gram of this model for unit filling fraction was studied in 
Ref. 52. 

It can be seen from the figure that incorporating the 
capacity for anyons to braid around one another results 
in an increase in the ground state energy per particle. 
This fact is reminiscent of the property that a system of 
identical fermions have a higher energy than bosons due 
to Pauli exclusion principle in real space. This also im¬ 
plies that there might exist a Pauli-like exclusion princi¬ 
ple for anyons too, at least in some regimes.®^ We also see 
from the figure that while the bosons and fermions have 
a paricle-hole symmetry which is reflected in the sym¬ 
metric ground state energy around half-filling v = ^, the 
system of Fibonacci anyons on the ladder does not dis¬ 
play this symmetry. One of the consequences of particle- 
hole symmetry is that the ground state energies E at 
filling fractions v and \ — v should be equal. While this 
is known for fermions and bosons, and reproduced by our 
numerical results as shown in Table. H, we see from our 


V 

Ahcb 

FnCFib 

Esf 

0 

0 

0 

0 

1/8 

-0.71162 

-0.70397 

-0.70015 

2/8 

-1.26597 

-1.19102 

-1.13658 

3/8 

-1.61707 

-1.44620 

-1.35273 

4/8 

-1.74300 

-1.52085 

-1.43534 

5/8 

-1.61707 

-1.41803 

-1.35271 

6/8 

-1.26597 

-1.15486 

-1.13660 

7/8 

-0.71162 

-0.68857 

-0.70015 

1 

0 

0 

0 


TABLE II. The values of the ground state energy E at var¬ 
ious filling fractions v corresponding to figure Fig. 9. The 
subscripts in i?(,) are “HCB” for hardcore bosons, “HCFib” 
for hardcore Fibonacci anyons, and “SF” for spinless fermions. 
The values are given to five decimal places. The ground state 
energies of bosons and fermions are symmetric around half¬ 
filling, but not so for Fibonacci anyons. 

numerical results that this no longer holds for some non- 
abelian anyon model such as Fibonacci anyons, though 
in this instance the breakdown of particle-hole duality is 
weak in the sense that it has only a very small impact on 
ground state energies. Interference of braiding particles 
raises the ground state energies, and thus the higher fill¬ 
ing fractions v> 1/2, e.g. v = 5/8, have slightly higher 
energies than the \ — v states, e.g. v = 3/8. 

The origin of the breakdown in the particle-hole dual¬ 
ity is in the difference in the fusion degrees of freedom 
of the particle types. For systems of bosons or fermions, 
the fusion space is one-dimensional, independent of the 
number of particles. In contrast, for non-Abelian mod¬ 
els such as the Fibonacci model, the fusion space grows 
exponentially with the number of anyons and hence is 
not symmetric under particle hole exchange. Braiding 
acts non-trivially on the fusion degrees of freedom and 
changes the ground state energy in way that is not par¬ 
ticle hole symmetric. 

2. Phase diagram of the Golden Ladder 

We further test our ansatz by studying the entangle¬ 
ment structures of ground states of interacting Fibonacci 
anyons on the ladder at unit filling. This model has been 
studied in Ref. 52, and we verify our ansatz by reproduc¬ 
ing known phases of the model at specific values of the 
tunable parameters. At unit filling, there is a single lo¬ 
calized Fibonacci anyons per site of the ladder and there¬ 
fore hopping rates are everywhere zero. This is a quasi- 
ID generalisation of the Golden Chain,which might 
be called the Golden Ladder. The relative interaction 
strengths of the legs and rungs of the ladder, including 
both antiferromagnetic and ferromagnetic couplings, may 
be parameterized on a circle (see Fig. 10, where the ferro¬ 
magnetic or antiferromagnetic natures of the interactions 
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FIG. 10. The horizontal and vertical interaction stengths 
(JhjJv) on the legs and rungs of the ladder are parameter¬ 
ized by (cos 6, sin 9) along the legs and rungs respectively. The 
labels within each quadrant indicate the nature of the inter¬ 
action, whether antiferromagnetic or ferromagnetic. 



FIG. 11. Scaling of entanglement entropy 5 as a function 
of block size r, for different angles on the circle (6), which 
correspond to different ratios of coupling strength between 
the legs and rungs of the ladder. 


in each sector are indicated). 

We evolve this model to ground state using TEBD, and 
compute the scaling of the block entanglement entropy 
from von Neumann’s relation, 

S{r) = -Ti{pr log pr), (32) 


where pr is in general the reduced density matrix of a 
block of r sites, here r rungs. From conformal field the¬ 
ory, the scaling of entanglement entropy on a system with 
an open boundary is 


S'(r) = |logr, (33) 

where c is the central charge of the system at criticality. 
This relation means that, for a critical model, the entan¬ 
glement block scaling—computed from the MPS ground 
state representation—should display a logarithmic rela¬ 
tion with the block size. The central charge c can then 
be extracted from the relationship 


3 -^(^ 2 ) - S{rx) 
log t -2 - log Ti ■ 


(34) 


The block entanglement entropy for various parameter 
regimes are shown in Fig. II, and their central charges 
are indicated in Fig. 12. As seen in Fig. II the finite 
bond dimension of the MPS causes entanglement to ar¬ 
tificially plateau over larger distances r = 1^2—7-11, but 
calculation of c using Fq. (34) may be performed for any 
separation r prior to this plateau, where an appropriate 
linear correlation is obtained between S{r 2 ) — S{ri) and 
log r 2 - log ri. 

One can interpret this Fig. 12 by considering how the 
physics of the interacting Fibonacci anyon changes as the 
parameterization angle 6 is varied. When 0 = 0, there 
are no couplings along the rungs and we have 2 chains 
of Fibonacci anyons with antiferromagnetic interactions. 
The system in this parameter regime is gapless and has 
a central charge which is twice that of a single chain. 



FIG. 12. The central charge of the underlying CFT extracted 
from the scaling of the entanglement entropy of Fig. 11 are 
shown at the parameter points we considered. When ver¬ 
tical coupling is set to zero and A range from —1 to -|-1, 
we obtain central charge which doubles that of single critical 
FM or AFM chains which lies on the equator. The values of 
cafm Chain and Cfm Chain are given as cafm chain = 0.704 and 
cfm Chain = 0.801. Phase boundaries for this model may be 
found in Fig. 11(a) of Ref. 52. 


i.e. 2 X 7/10. Fven though the MPS most naturally 
yields exponentially decaying correlators, we are never¬ 
theless able to extract an approximate value for the cen¬ 
tral charge, c = 1.405, from the linear part of the curve. 
When 0 = 7r/4 and 0 = 37r/4, the vertical couplings are 
antiferromagnetic favouring pairs of Fibonacci anyon fus¬ 
ing into the vacuum charge. This phase is gapped with 
central charge c = 0. When 0 = 7r/2, the Hamiltonian 
favours fusion of pair of r charges on each rung to the vac¬ 
uum charge, and is hence a product state which is unique 
and gapped. The phase is not critical and has a central 
charge c = 0. At the 0 = tt point, the horizontal coupling 
A = — 1 is ferromagnetic, while the vertical coupling A 
is zero, and the ladder reduces to two copies of a ferro¬ 
magnetic Golden Chain. From our numerical simulation. 
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e 

Csim 

CTheo. 

0 

1.405 

2 X 7/10 

7 r /4 

0 

0 

'K j'l 

0 

0 

37 r /4 

0 

0 

TT 

1.629 

2 X 4/5 

j A: 

0.801 

4/5 

Z-kII 

0 

0 


0.704 

7/10 


TABLE III. The table shows the obtained central charges 
csim. from our numerical simulations, compared against their 
theoretical values cxheo. known from conformal field theory, 
at the interaction strengths parameterized by 6 according to 
Fig. 10. The values correspond to the points shown in Fig. 12. 
Where a model is not critical, and hence not described by 
CFT, we have substituted their central charge c with zero. 

we computed a central charge of c = 1.629 which is close 
to the expected theoretical value of c = 2 x 4/5. At the 
point 9 = 57r/4, the horizontal and vertical couplings are 
ferromagnetic. Fusion of the r charges on the rungs and 
legs favours projection to the r channel (triplet state). 
This can easily be pictured by considering a linearized 
version of the ladder. Nearest neighbour r charges on 
the rung becomes nearest neighbour on the chain and 
nearest neighbours on the legs becomes next-to-nearest 
neighbour on the chain.Heuristically, fusion to the r 
fusion channel makes the ladder effectively like a single 
Fibonacci chain and therefore has the same central charge 
as a single chain. We obtain a central charge of c = 0.801 
which is close to the expected c = 4/5. When 6 = the 
vertical coupling J„ = — 1 is ferromagnetic while the hor¬ 
izontal coupling is zero. This favours projection of neigh¬ 
bouring r charges on the rungs into the r channel. The 
ladder reduces to a chain of decoupled r charges which 
has an exponentially large degeneracy in intermediate fu¬ 
sion degrees of freedom. Hence a generic ground state 
at this point obeys a volume law rather than area law. 
This system is gapped and not described by conformal 
field theory. The parameter point 9 = lir/A correspond 
to horizontal antiferromagnetic coupling on the leg and 
vertical ferromagnetic coupling which is effectively an an¬ 


tiferromagnetic interacting chain. The obtained central 
charge is c = 0.704, being close to the expected value of 
c = 7/10. Our findings are in agreement with known re¬ 
sults showing that the entire upper semicircle is gapped 
while the lower semicircle is gapless with the exception 
of the indicated point at 9 = 37r/2. 

Table III compares the extracted central charges with 
their expected theoretical values. 

VI. CONCLUSION 

In this paper we show how the anyonic tensor net¬ 
work formalism of Refs. 40, 43, and 44 may be ap¬ 
plied in the context of particles admitting multiple 
charge labels, specifically an anyonic charge and a U(l) 
charge, here corresponding to particle number. We con¬ 
structed test models involving both hopping and interac¬ 
tion terms, with this construction being explicitly elabo¬ 
rated in the Appendix. Application of the anyonic in¬ 
finite TEBD algorithm'^^ permitted calculation of the 
ground states of these systems, their entanglement en¬ 
tropies, and central charges. In doing so, we successfully 
reproduced elements of the phase diagrams for these sys¬ 
tems which have previously been obtained using exact 
diagonalisation.^^’^^’^®^^^ 

This paper consequently demonstrates the feasibility of 
applying anyonic TEBD to systems of particles admitting 
both anyonic and U(l) conserved charges. The method 
presented here can be used to probe new regimes of the 
physics of anyons such as equilibrium phases of quasi-ID 
systems of braiding anyons at arbitrary density as well as 
non-equilibrium dynamics of anyons in two dimensional 
systems at low density. The later could be used to study 
the robustness of large size topological quantum com¬ 
puters/memories to errors induced by coherent propaga¬ 
tion of erroneous anyons created by thermal fluctuations 
which braid around logical degrees of freedom. 
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A system of anyons which does not fuse to the total vacuum 
charge I can be considered to be within a larger system 
with a nontrivial boundary charge that enforces a complete 
annihilation of the anyonic charges to the vacuum charge. 


Appendix A: Basics of Anyonic Tensor Networks 

Using tensor networks to simulate quantum systems 
involves choosing a network pattern of connected tensors 
along with a choice of an algorithm that optimizes the 
representation of the many body state.In what follows, 
we review the basic objects common to most anyonic ten¬ 
sor networks (TN). In the main text, we used these ob¬ 
jects to construct the anyonic MPS, and our modified 
Anyonic-U(l) MPS ansatz. 


1. Components of anyonic tensor networks 

In the discussion below, we assume some familiarity 
with theory of anyons as described in Refs. 49, 55, and 56. 
For a pedagogical introduction to anyons, see e.g. Refs. 
21 and 57. 

The basic objects in any tensor network include vectors 
(or one-index tensors), matrices (or two-index tensors), 
and more generally, n-index tensors. We start by exam¬ 
ining how the basis of states for a system of anyons can 
be enumerated, and how it is used in constructing the 
anyonic equivalents of the above-named TN objects. 
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a. Anyonic basis enumeration and fusion lookup tables 

The basis of the Hilbert space of anyons is described by 
a labeled directed fusion tree (see Fig. 1) where the charge 
c on any incoming edge at a vertex is determined from 
the charges a and b of the two outgoing edges around the 
same vertex, according to the fusion rules of the anyon 
theory 

(Al) 


which implies that charges a and b are allowed to fuse 
to possibly several different charges c. For example, in 
the fusion tree in Fig. l(ii) charges bi and 03 fuse into all 
possible charges 62 • We restrict to multiplicity-free anyon 
models in this paper, i.e. G {0,1}, which includes 
some of the models most relevant to current experiment 
such as Ising anyons and Fibonacci anyons. The itera¬ 
tive fusion process described by a fusion tree makes the 
mathematical description of the collective state space to 
be non-local and hence not naturally reducible to a tensor 
product space of individual anyonic degrees of freedom. 

The description of labeled fusion trees can become ex¬ 
tremely verbose in the limit of large system sizes. The 
amount of data needed to specify the labeling of a fusion 
tree can be greatly reduced if one only enumerates the la¬ 
belings having a particular total charge at the trunk. To 
this end, let c be total charge at the trunk of the fusion 
tree [see Fig. l(i)] and assign an index /ic = 1, 2, • • • , z/c 
in increasing numerical order to each unique labeling of 
the fusion tree. The index pLc is called the degeneracy in¬ 
dex and r'c is called the degeneracy of charge All the 
fusion trees are therefore concisely labeled by the multi¬ 
index 7 = (c, /ic), with c as the total charge label and fj-c 
as its degeneracy index. 

This assignment can be described using a “fusion 
lookup table” where the sets of charges corresponding to 
unique labelings of the fusion tree are recorded in rows, 
and each row is additionally labeled with (i) the total 
charge c and (ii) an index value /Tc such that each pair 
(c, /ic) is unique. The use of fusion lookup tables serve a 
two-fold purpose in anyonic tensor networks. First they 
can serve as a “cache”, allowing for a total recovery of 
all the charges labeling a particular fusion tree which 
may be useful when performing diagrammatic (topologi¬ 
cal) manipulations of the fusion tree. Secondly, they can 
be used to significantly reduce the computational cost of 
contracting an anyonic TN, a fact which has previously 
been discussed in the context of non-anyonic symmetric 
TNs (see Refs. 43, 44, and 46 for more details). When 
constructing a fusion table, it is required that there be 
a one-to-one correspondence between the charge label¬ 
ings of the fusion tree and the charge/ degeneracy multi¬ 
indices which are assigned to these labelings. 




FIG. 13. (i) Anyonic state vector and (ii) its Hermitian con¬ 
jugate. 


b. Anyonic state vector 

An anyonic quantum state jlk) can be written as a 
weighted superposition of all labelings of a fusion tree 
having a total vacuum charge I.®® More compactly, in the 
multi-index notation, the quantum state can be written 
as 

(A2) 

7 

where 7 = (I, /ii) is an index enumerating all the valid 
fusion trees. If all the enumerated fusion trees are associ¬ 
ated with normalized anyonic diagrams, which is referred 
to as the implicit normalization scheme, then the state 
amplitudes can be arranged as a column vector in 
the standard basis. Following the diagrammatic nota¬ 
tions employed for anyonic tensors in Ref. 40, we depict 
the anyonic quantum state by a filled circle with a central 
leg enumerating all the multi-indexed bases, and an unla¬ 
beled tree structure, as shown in Fig. 13(i). It should be 
noted that if topological manipulations were to be per¬ 
formed on the fusion tree, such as vertically bending a 
line opposite to its orientation, it is preferable to use dia¬ 
grammatic isotopy convention. In such a case, we adopt 
the prescription given in Ref. 49, where the fusion dia¬ 
grams are weighted with certain pre-factors of quantum 
dimensions of the anyonic charges on the fusion tree. For 
the fusion basis shown in Fig. 1, the pre-factor would be 
/ . \i/4 

I 7— j ° j — . These normalization factors are then 

also absorbed into the amplitudes defining the anyonic 
state vector. This is referred to as the explicit normaliza¬ 
tion scheme. During topological manipulations, all the 
charges labeling a particular fusion tree can be recovered 
from the fusion lookup tables and used in computing the 
necessary data associated to that operation. We work 
exclusively in the explicit normalization scheme, where 
each vertex is normalized according to diagrammatic iso¬ 
topic convention. The differences between working in 
implicit and explicit normalization scheme, collectively 
called mixed normalization, are treated in the recent any¬ 
onic DMRG paper. 
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The Hermitian conjugate of the state written as 

(A3) 

7 

is represented diagrammatically as in Fig. 13(ii) where 
the unlabeled fusion tree is reflected vertically and all its 
arrows are reversed. The coefficients of the vector are 
also complex-conjugated. 

c. Anyonic Matrix Operator 

In conventional quantum theory, an operator O : V —> 
Y' is written in the bra-ket notation as, 

0 = Y.OfAj'){3\- (A4) 

j'J 

where the indices j and j' enumerates basis states in V 
and Y' . Example of such operators include Hamiltonians, 
density matrices, projectors, etc. 

In a similar vein, an anyonic operator acting on a set 
of anyonic charges with total charge c does not change 
the total charge. The operator Oc ■ ^ 

takes states of anyons oi, 02 , • • • , Ofc to states 
of anyons , 02 , • • • , aj,/ without changing the conserved 
total charge c. As such, the operator O = 0^ Oc can be 
constructed as a block-diagonal matrix with each block 
indexed by the conserved anyonic charge c. Each block 
matrix Oc is constructed by enumerating (as in Fig. 1) all 
the fusion tree bases fusing to that charge. As such the 
charge-conserving matrix is indexed by the multi-index 
7 = {c, yLc) for fusion trees and 7 ' = (c,pl'c) for splitting 
trees. The anyonic operator can therefore be written as 

Oc = E^7',7 IV) (7h (A5) 

7'O' 

where 7 = (c, jic) and 7 ' = (c, implying charge con¬ 
servation. The matrix elements will depend on the partic¬ 
ular physics of the system. The anyonic matrix operator 
is represented diagrammatically by Fig. I4(i), where the 
midti-indices 7 = (c, /ic) and 7 ' = (c, /i(,) enumerate all 
the fusion and splitting trees. The vertex normalization 
factors of the fusion/splitting trees are absorbed into the 
matrix operator. 

d. Anyonic rank-3 tensor 

The anyonic matrix operator can be extended to a 
rank-3 tensor where the tensor elements are indexed 
by three multi-indices a, /3 and 7 . An anyonic rank-3 
tensor T“’^ may be represented in the manner shown 
in Fig. I4(ii) where the leaves on each branch of the 
tree are enumerated and assigned a multi-index nota¬ 
tion a = {a,tLa), j 3 = {b, pib) and 7 = (c,/ic)- All the 



FIG. 14. (i) The anyonic matrix operator O and (ii) the 

anyonic rank-3 tensor T. The tree has been normalization 
according to diagrammatic isotopy. 

vertices on the leaves fulfill the fusion rules during enu¬ 
meration of the basis, and the implicit vertex contained 
within the grey circle also obeys the fusion rules of the 
anyon model. The tensor T“’^ is indexed by 7 = (c,/ic), 
a = {a, tia) and /3 = {b^fib), where the charge triplet 
(a, 6 , c), obtained from each of the subtrees, has to be 
compatible with the orientation of the vertex on the ten¬ 
sor, i.e. a X b ^ c. The explicit form of the tensor T 
is 

Va,Vb,Vc 

T = ^Kb E {T:X:’''"\k^akib){k^c\. (A 6 ) 

a,b,C fla,nb,nc = t 

The direct sum implies that tensor T is composed block- 
wise from tensors indexed by the charges of the subtrees, 
with each block then being indexed by the degeneracy 
index of the compatible fusion trees. 

There are more objects that can be implemented to 
manipulate anyonic tensor networks^*^’^^’^^, but as the 
MPS is a trivalent tensor network, the anyonic tensors 
we have reviewed are sufficient to construct the anyonic 
MPS. 

Appendix B: Anyon Model Data 

An anyon model is minimally specified by the follow¬ 
ing data: a set of charges A, fusion rules for the charges 
the braid matrix i?, and the F-tensor F. All other 
quantities can be derived from these data. The anyon 
models we used for testing our numerical method in this 
paper are Fibonacci anyons, Z 2 (spinless) fermions and 
Zoo bosons, but the charge spectrum of the physical site 
is restricted to Ahcb = { 0 , 1 } corresponding to the hard¬ 
core constraint. 


e. Fibonacci anyon data 

The Fibonacci anyon model consists of two charges, 
vacuum (I) and Fibonacci anyon (r). Hence A = (I, rj 
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where the charges have quantum dimensions, di = 1, 
dr = respectively. The fusion rules obeyed by 

charges are 

1x1 = I, IXT = TXl = T, TXT=I + T. (Bl) 

The fusion tensor N has components = 0 when a x 
b ^ c for all a,b,c G A. The nonzero components are 
given by 

iv/i = iv;i = iVi; = ivi, = iv;, = i. (b2) 

The i?-inatrix has nonzero components 

^ g-4,r*/5^ j^rr ^ ^ pAr ^ prl = i?n = 

(B3) 

for compatible charges and zero otherwise. The nontriv¬ 
ial f-move coefficients are 



where (j) = , and e,/ S {I, t}. The remaining F- 

move coefficients are given by 

(Ff= N:,NlNtNtf. (B5) 

/. Fermions and Bosons Data 

Fermions and bosons can be studied within the the¬ 
ory of anyons and consequently using anyonic tensor net¬ 
works such as the anyonic MBS. The wave functions of 
fermions and bosons acquire phase factors of —1 and -1-1 
respectively under particle pair exchange. 

The particle spectrum of bosons is the set of pos¬ 
itive integer charges, denoted as Zoo = {0, 1 , 2 ,.. .} with 
the fusion rule being ordinary addition, while the charge 
spectrum of fermions is Z 2 = {0,1} with fusion rules cor¬ 
responding to addition modulo 2, so 1 x 1 —>■ 0. The 
quantum dimensions are trivial with do = di = 1 for 


fermions and d^ = 1 V g for bosons. The fermionic ex¬ 
change factors (permutation factors) are encoded in the 
fermionic i?-matrix as 

i?oo = Rl° = = 1, Rl^ = -1, (B6) 

while for bosons all valid entries are trivially equal to 
one. The F-matrix for both particle types fulfills 

(^ahc)/ ^ N:,NiXcN^f (B7) 

Appendix C: Derivation of the anyonic t-J 
Hamiltonian on a ladder 

We now give the explicit derivation of the Hamiltonian 
d on a plaquette. This Hamiltonian consists of itinerant 
and (Heisenberg) interaction terms written as 

h = —tfi (ha-^c=I + bi,^ii=I + h-C^ -|- Jh -|- Hj 

~ + h-C^ -I- — -I- . 

(Cl) 

We derive the matrix representation of each term of the 
Hamiltonian in the first basis {|/)} shown in Fig. 8(i). 
For operators like 6 q_>c=i, H^^, etc., which couple anyons 
on the legs of the plaquette, we transform the basis |/) 
to the basis \II) in Fig. 8(ii), derive the action of the 
Hamiltonian in the basis \II), then transform back to 
the basis |/). 

Therefore we first show how the two bases transform, 
and later show the derivation of the Hamiltonian for the 
plaquette. 

1. Fusion Tree Basis Transformation 

The transformation between the two chosen bases of 
Fig. 8 are obtained as follows: 


abed 

7 


abed 





a c b d 


a 


V 

7 

,-n'' 


= E,,. {FrT' FF)! ^5° W 


a c b d 


7 




, -1 



7 


h d 





(C2) 
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The above equation can be written more succinctly as 
abed a c b d 



a 





1 7 

where the tensor is defined according to 


(C3) 


K,X 


X (^ 7 “), ■ 

(C4) 

Using Dirac bra-ket notation, Eq. C2 can alternatively 
be written as 

\{ab;a){cd;/3){a/3-,^)) = ^ |(ac; K)( 6 d; A)(kA; 7 )). 

K,X 

(C5) 


2. Anyonic t-J Hamiltonian on a plaquette 

The anyonic local Hamiltonian h on a plaquette is 
given by 

h = -th {ha^c=i + bb^d=i + h-c) - Jh (n^,c + 

- I (b^^b=i + k^d=i + h-c) - . 

(C 6 ) 

Whereas charge label a may take any value from the par¬ 
ticle spectrum, the vacuum charge will be denoted by I 
and a nontrivial anyonic charge by qq. For example, in 
the Fibonacci anyon theory, oq = t. The derivation is 
quite general and can be used with any anyon model. 
Note that numerical factors such as vertex normalization 
factors and loop factors are not accounted for here. We 
account for them during the implementation of the any¬ 
onic TEBD algorithm. 

We proceed by first deriving all the kinetic energy 
terms and then derive all the interaction terms similarly. 
All the operators in the Hamiltonian are applied to the 
fusion tree on the left hand of Eq. C2 which is represented 
in Dirac notation in Eq. C5. 

Kinetic terms : The terms contributing to the kinetic 
energy are the braid operators, whose matrix elements 
are derived below. 

(i) The matrix element of the braid operator ba->.b=i is 
given by 

(^a—l) ^a^ao^b,!: ^a\b^b',a^c'^c^d\d^a'p',P ^ 

where we have used, for the sake of conciseness, the no¬ 
tation {ba^b=i) as a shorthand for 

{{a'b'; a'){c'd'\'y)ba^b=iiab] a){cd-, (3)(a(3; 7 )). 

(C 8 ) 


(ii) The matrix element of the braid operator bc^d=i 
is given by 

{bc^d—l} ^c,co^d,I: ^a\a^b',b^c',d^d'^c^a\a^p',p- 

(iii) The matrix element of the operator ba^c=i in¬ 
volves braiding of anyonic charge a with b. The charge 
c has to be vacuum for the process to have a nonzero 
amplitude. Its action on the basis |(a&; a)(cd;/3)(a/3; 7 )) 
is given by 


ba^c=i |(a 6 ;a)(cd;/ 3 )(a/ 3 ; 7 )) 

= Ek,a Qapk^c=i\{ac; K){bd- A)(kA; 7 )), 
= J2i^,\Qap^a,aoSc,i\{ca; K){bd; A)(kA;7)) . 

The expectation value ( 6 a_>.c=i) is 


{ba^c=i) = {ia'c';K'){b'd'-,X'){K'X';-f)\Q*j'pk 

K,X 

^^a,ao^c,I \{ca]K){bd] A)(kA; 7 )), 

(CIO) 

which simplifies to 


(ba^c^l) =J2Q':,^piQk"x'Sa,aoSc,lSa',cSc',aSb',bSd',d. 

K,X 

(Cll) 

(iv) The expectation value {bb^d=i) is similarly given 
by 

(bb^d=l) = '^^Q'kp{Qk^\ ^b,boSd,l5a' ,abc' ,cSb' ,dbd' ,b- 

k,A 

(C12) 

Interaction terms : The interaction terms consist of 
projectors whose matrix elements are derived similarly 
to the braid terms. The projection favours fusion of non¬ 
trivial anyons to the vacuum charge. 

i) The action of the projector H^ ^ on the fusion basis 
is given as 

nL |(a 6 ;a)(cd;/ 3 )(a/ 3 ; 7 )) = |(a 6 ; a)(c(I;/3)(a/3; 7 )) 

(C13) 

where the element n“j = 1 if a = I (vacuum) and a = ag, 
b = bg (i.e. nontrivial charges). The expectation value of 
the projector 11 *^ is thus 


(lifjf, ) — ba\ao^b\bo^a\abb\bba' ,cbd\db Id ,P* ■ 

(CM) 

(ii) The matrix element H),^ of the projector is similarly 
given as 

) — ba\abb' ,bba' ,a^c,CQb(l^fiQb0^gSc\c^d' ,d^P,P' • 

(C15) 

(iii) The action of the projector HJ)^ on the basis 
|(a6;Q;)(cd;/3)(a/3;7)) is 


nL \iab-,a){cd]P){al3;j)) 

= EQS^n:j(ac;«)(6d;A)(«A;7)) 

K,X 

= Ek,a Qap^a,aoSc,coS,^,i l(ac; K){bd; X)(kX; 7 )). 
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The matrix element (IlJjc) is 

\ 


c,cq ^a' ,a^c 


c^b'b^d' 


d- 


(C16) 


(iv) The matrix element (flj^) is similarly given by 

K,X 

(C17) 


